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RATIONAL  PULSE-SHAPING  FILTERS  FOR 
ALMOST  MSK-LIKE  MODULATION 


1.  INTRODUCTION 

MSK-like  modulation  formats  are  those  for  which  (a)  the  instantaneous  phase  of  the  transmitted 
signal  changes  by  exactly  90°  in  each  bit  interval,  and  (b)  the  instantaneous  amplitude  (or  envelope)  of 
the  transmitted  signal  is  constant.  Well-known  members  of  this  family  include  MSK  [1]  (minimum- 
shift  keying),  SFSK  (sinusoidal  frequency-shift  keying)  [2],  and  DSFSK  (doubly  sinusoidal  frequency- 
shift  keying)  [3],  These  formats  are  used  where  a  constant  envelope  is  needed  but  with  a  sharper  spec¬ 
tral  rolloff  than  is  provided  by  rectangular  QPSK  (quadriphase-shift  keying).  Rectangular  QPSK  rolls 
off  asymptotically  at  6  dB/octave,  MSK  rolls  off  at  12  dB/ octave,  SFSK  at  24  dB/octave,  and  DSFSK  at 
36  dB/octave.  The  width  of  the  main  spectral  lobe  is  greater  for  the  formats  with  the  sharper  asymp¬ 
totic  spectral  rolloffs.  Of  these  formats,  only  MSK  is  easily  generated,  and  consequently  MSK  is  used 
regularly  both  for  data  communication  and  for  radar  pulse  compression  with  binary-coded  sequences. 
The  other  formats,  desirable  for  their  lower  spectral  sidelobes,  generally  require  complex  filters  or  pre¬ 
cision  voltage-controlled  oscillators,  making  them  difficult  to  implement,  or  rely  on  digital  signal- 
generation  methods,  making  them  unsuitable  for  high-speed  operation. 

In  this  report,  simple  implementations  are  sought  not  by  looking  for  simple  approximations  to 
existing  modulation  formats  in  the  MSK-like  class,  but  by  searching  within  a  particular  easily  imple¬ 
mented  class  for  those  modulation  formats  that  come  closest  to  being  MSK-like.  This  approach  resulted 
in  several  new  modulation  formats  that  (a)  change  phase  by  90°  in  each  bit  interval;  (b)  have  almost 
constant  envelope;  and  (c)  can  be  implemented  to  good  approximation  using  low-order  rational  transfer 
functions  to  filter  in-phase  and  quadrature  NRZ  (nonreturn-to-zero)  waveforms.  These  new  formats 
were  derived  by  First  starting  with  the  class  of  waveforms  that  can  be  generated  by  passing  in-phase  and 
quadrature  NRZ  waveforms  through  identical  filters  with  rational  transfer  functions,  modulating  the 
two  waveforms  onto  quadrature  carriers,  and  summing.  Further  restricting  the  class  of  waveforms  con¬ 
sidered  to  those  that  change  phase  by  90°  in  each  bit  interval,  the  locations  of  the  singularities  of  the 
transfer  functions  of  the  pulse-shaping  filters  were  optimized  with  respect  to  the  mean  of  the  fourth 
power  of  the  signal  amplitude  with  signal  energy  constrained  to  be  constant.  This  resulted  in  modula¬ 
tion  formats  with  very  small  variations  in  amplitude. 

If  small  amplitude  fluctuations  were  not  present,  these  modulation  formats  would  be  MSK-like.  I 
use  the  term  almost  MSK-like  to  refer  to  modulation  formats,  such  as  those  presented  here,  that  have 
such  small  amplitude  variations  that  they  can  be  treated  as  MSK-like  for  all  practical  purposes.  I  refer 
to  the  optimum  almost-MSK-like  modulation  formats  presented  in  this  report  as  RSK,  for  rational-shift 
keying.  The  term  rational  emphasizes  that  these  formats,  unlike  the  MSK-like  formats,  can  be  gen¬ 
erated  with  rational  pulse-shaping  filters.  The  -shift  keying  emphasizes  the  similarity  to  the  MSK-like 
formats.  The  term  is  more  mnemonic  than  accurate,  however,  for  the  instantaneous  frequency  shift  of 
RSK  signals  cannot  be  described  as  a  rational  function.  Nevertheless,  I  use  the  term  RSK  here  for  lack 
of  a  more  descriptive,  equally  concise  term.  To  differentiate  between  the  various  RSK  formats 
described  in  this  report  I  will  use  the  order  of  the  associated  pulse-shaping  filter.  This  report  derives 
fourth-,  sixth-,  and  eighth-order  RSK  formats. 


Manuscript  approved  November  5,  1984. 
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The  remainder  of  this  report  is  organized  into  five  sections.  The  section  following  this  introduc¬ 
tion  establishes  notation  and  provides  the  minimum  background  necessary  to  understand  the  remainder 
of  the  report.  The  class  of  modulation  formats  optimized  is  then  developed.  The  optimization  criterion 
used  to  choose  the  best  formats  within  that  class  is  then  presented,  followed  by  a  presentation  of  the 
RSK  formats  that  resulted  from  the  optimization  process.  The  summary  includes  the  RSK  parameters 
in  tabular  format  for  convenient  reference.  Various  lengthy  mathematical  details  are  relegated  to 
appendixes. 


2.  BACKGROUND:  QUADRATURE  DATA-MODULATION  FORMATS 

In  this  section  I  describe  several  related  modulation  types  as  a  vehicle  for  establishing  notation.  I 
first  describe  a  very  general  modulation  format,  and  then  I  describe  various  subsets  of  that  format  until 
I  have  described  the  class  of  modulation  formats  optimized  in  the  following  sections. 

QAM 

Consider  amplitude  modulating  a  carrier  with  a  complex  modulation  function  mit),  where* 

mit)  -  iit)  +  jq(t), 

with  iit )  and  q{t )  real,  and 

sc  ( t )  =  Re{  m  ( t )  eJo>c  r}  =  iit)  cos  wcf  —  qit)  sinwcf.  (1) 

This  is  referred  to  as  quadrature  amplitude  modulation  or  QAM.  The  complex  modulation  function  m{t) 
can  be  recovered  from  scit)  if  mit )  is  band-limited  to  frequencies  co  >  —  <wc.  Conversely,  any  real 
scit)  with  no  DC  component  can  be  generated  with  a  suitable  m(t).  Thus,  QAM  is  a  very  general 
modulation  format;  all  of  the  modulation  formats  discussed  in  this  report  are  special  forms  of  QAM. 

QPSK 

QPSK,  or  quadriphase-shift  keying,  is  a  form  of  QAM  in  which 

mit )  =  AeJ<t>  £ amp(t  —  2 mT),  (2) 

m 

with  each  am€  { 1 1  ,—j)  and  with  pit),  termed  the  baseband  pulse  Junction,  satisfying 

T)p(d  +  2kT)  =  S*  '  (3) 

(8k  is  the  Kronecker  delta  function)  for  some  fixed  delay  d  and  scale  factor  tj.  The  latter  condition  is 
the  time-domain  equivalent  of  Nyquist’s  second  criterion  and  is  required  to  prevent  Intersymbol 
interference,  i.e.,  to  ensure  that  at  times  t  =  d  +  2kT  only  one  term  contributes  to  the  sum  in  Eq.  (2). 
If  pit)  is  a  rectangular  pulse  of  width  2T,  so-called  rectangular  QPSK  results.  The  term  QPSK  itself  is 
reserved  by  some  authors  for  rectangular  QPSK. 

QPSK  is  frequently  used  in  data  communication,  where  bits  are  made  available  for  transmission  at 
rate  T-1.  Bits  are  grouped  into  pairs,  and  each  pair  is  mapped  into  one  of  the  am,  with  each  of  the  four 
permutations  of  two  bits  mapped  into  one  of  the  four  possible  values  for  the  am.  Rectangular  QPSK  is 
also  occasionally  used  in  radar  pulse  compression,  where  the  quaternary  sequence  {«m}  is  carefully 
chosen  to  have  desirable  autocorrelation  properties.  There  is  no  reason  why  nonrectangular  QPSK 
could  not  be  used  in  pulse  compression  as  long  as  the  amplitude  of  the  baseband  pulse  function  pit) 
decayed  quickly  in  either  direction  from  the  peak. 


’The  unit  imaginary  quantity  is  represented  in  this  report  by  j  rather  than  /.  This  is  common  in  the  literature  of  electrical 
engineering  and  here  permits  the  conventional  use  of  /  for  in-phase. 
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The  conventional  implementation  of  QPSK  effectively  rearranges  Eq.  (2)  to  obtain 


mit)  = 


J^amg(t  -  2 mT) 


®fit). 


(4) 


where  ©  indicates  convolution,  and  where  git )  ®fit)  =  pit)  of  Eq.  (2).  The  gate  function  git)  is 
defined  by 

git)  =  Uit)  -  Uit  -  2T),  (5) 


where  Uit)  is  the  unit-step  function,  i.e.,  git)  is  equal  to  unity  for  0  <  t  <  2T  and  zero  everywhere 
else.  Thus,  the  bracketed  sum  in  Eq.  (4)  is  just  a  pair  (real  and  imaginary  parts)  of  NRZ  digital 
waveforms.  Since  fit)  is  nearly  always  a  purely  real  function,  the  QPSK  implementation  correspond¬ 
ing  to  this  formulation  is  as  shown  in  Fig.  1.  The  convolution  of  Eq.  (4)  is  implemented  by  passing  the 
two  NRZ  data  streams  derived  from  the  «m  through  pulse-shaping  filters  with  impulse  responses  fit). 
The  upper  and  lower  signal  paths  in  Fig.  1  are  referred  to  respectively  as  either  the  i  and  q  channels  or 
arms.  The  pulse-shaping  filters  in  the  /  and  q  channels  are  referred  to  in  the  remainder  of  this  report 
by  their  conventional  names,  the  i  and  q  arm  filters. 


i  DATA 


I  2T  j 

I  I 


1 


I 

q  DATA 

'  i 


—  sin  wct 

Fig.  1  —  Conventional  QPSK  implementation 


OQPSK 


Attempts  in  the  data-communication  community  to  develop  modulation  formats  less  susceptible 
to  degradation  in  nonlinear  transmitters  led  to  the  development  of  staggered  or  offset  QPSK,  abbreviated 
OQPSK.  OQPSK  is  identical  to  QPSK  except  that  the  imaginary  part  q(t)  of  the  modulation  function 
m(t)  in  Eq.  (2)  has  been  offset  in  time  from  the  real  part  /(/)  by  time  T.  This  reduces  the  extreme 
envelope  fluctuations  in  sc(t)  of  Eq.  (1)  suffered  by  ordinary  QPSK  (especially  if  nearly  rectangular) 
whenever  am  =  -am_i  in  Eq.  (2). 


A  mathematical  formulation  of  OQPSK  in  terms  of  the  m(t)  of  Eq.  (2)  is  clumsy.  OQPSK  is 
more  concisely  defined  as  QAM  with 


mit)  =  '£/3mjmpit  -  mT), 

m 


(6) 


where  each  /3m€  {1,-1 },  and  where  pit)  is  constrained  as  for  ordinary  QPSK  above.  The  {/3m)  in  Eq. 
(6)  might  typically  be  chosen  for  data  transmission  as 


Pm 


Pm- 1  to  send  a  1 
—  Pm-\  to  send  a  0, 


(7) 
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where  /30  is  chosen  arbitrarily.  OQPSK  is  not  any  less  suited  to  radar  pulse  compression  than  ordinary 
QPSK. 

The  factor  AeJ<i>  of  Eq.  (2)  was  omitted  in  Eq.  (6)  in  the  interest  of  brevity.  It  should  be  under¬ 
stood  that  OQPSK  and  any  of  the  modulation  formats  described  subsequently  can  have  an  arbitrary 
amplitude-and-phase  factor  added. 


The  OQPSK  formulation  of  Eq.  (6)  can  also  be  rewritten  to  explicitly  indicate  the  use  of  arm 
filters  on  NRZ  data  streams: 


m{t)  = 


I  (3mjmg(t-rnT) 


®/(0. 


(8) 


The  corresponding  implementation  is  given  in  Fig.  2.  The  only  difference  between  Figs.  2  and  1  is  that 
in  Fig.  2  the  q  data  stream  has  been  offset  from  the  /  data  stream  by  time  7\  The  data  values  fed  to 
the  arm  filters  in  Fig.  2  are  related  in  a  simple  way  to  the  {/3m}  of  Eq.  (8),  but  the  details  are 
unimportant  here. 


CD5  wei 


i  DATA 


-sin  wct 


Fig.  2  —  Conventional  offset-QPSK  implementation 


It  is  interesting  to  note  that  for  baseband  pulse  functions  satisfying 

f)p(d  +  mT)  =  8m  (9) 

for  some  fixed  delay  d  and  scale  factor  tj  (compare  to  the  Nyquist  criterion  of  Eq.  (3)), 
m(d  +  kT)  =  ±  jm(d  +  (k  —  1)D,  i.e.,  the  OQPSK  signal  undergoes  a  net  phase  change  of  exactly 
90°  in  each  bit  interval.  The  subset  of  OQPSK  discussed  in  the  next  section  will  satisfy  this  condition. 


MSK-Like  Formats 

What  considerations  are  involved  in  creating  a  constant-envelope  OQPSK  signal,  that  is,  in  mak¬ 
ing  |m(f)|2  =  1  in  Eq.  (6)  using  a  real  baseband  pulse  function  p(t)l  Answering  this  question  requires 
looking  at  the  properties  of  |m(r)|2.  Take  m(t )  as  given  by  Eq.  (6)  and  multiply  by  its  conjugate; 

M/)|2  =  m{t )  m*{t)  =  £  £j3j3  njm~npit  ~  mT)  pit  -  nT),  (10) 

m  n 

where  the  *  superscript  indicates  complex  conjugation.  Since  the  {f3m}  are  data-dependent,  and  thus 
unpredictable,  the  only  way  \m(t)\2  can  have  a  predictable  value  is  if  the  only  terms  that  contribute  to 
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the  double  sum  of  Eq.  (10)  are  those  for  which  m  =  n.  The  resulting  requirement  is  that 
pit  —  mT)  pit  —  nT)  —  0  for  all  t  and  m  ^  n.  With  some  simple  substitution  this  becomes 

pit )  pit  —  mT)  =  0;  m  ^  0, 

which  is  easily  seen  to  be  equivalent  to 

pit)  =  0;  U-  t0 1  >  T  (11) 

for  some  f0,  i.e.,  pit)  can  only  be  nonzero  over  an  interval  of  2 T  in  t.  Using  this  assumption  that  only 
terms  for  which  m  =  n  can  contribute  to  the  sum  of  Eq.  (10),  Eq.  (10)  becomes 

\mit)\2=  ^p2it-  mT).  (12) 

m 

Notice  that  this  is  periodic  with  period  T.  To  force  \mit)\2  to  unity  it  will,  therefore,  be  adequate  to 
force  it  to  unity  over  one  period.  Since  Eq.  (11)  is  assumed  satisfied,  only  two  terms  of  Eq.  (12)  con¬ 
tribute  at  any  one  time.  Focusing  in  on  one  particular  interval  of  width  T, 

\mit)\2  =  p2it)  +  p2it  +  T)\  t0-  T  <  t  <  t0.  (13) 

Equating  this  to  unity  produces  the  second  constant-envelope  requirement, 

p2it)  +  p2it  +70  —  1;  t0  -  T  <  t  <  t0.  (14) 

Equations  (11)  and  (14)  are  the  two  conditions  that  are  frequently  cited  as  defining  the  MSK-like  sub¬ 
set  of  the  OQPSK  modulation  formats  (often  with  t0  set  to  0).  The  most  well-known  modulation  for¬ 
mat  of  this  class  is  minimum-shift  keying,  or  MSK,  described  by  Eqs.  (6)  and  (7)  and 

pit)  =  cos  M  <  T. 

A  close  inspection  of  this  definition  reveals  that  MSK  is  also  a  special  case  of  FSK  (frequency-shift  key¬ 
ing),  with  scit)  of  Eq.  (1)  having  frequency  oc  +  tt /2T  while  a  data  1  is  being  transmitted  and  fre¬ 
quency  o)c  -  tt /2T  while  a  data  0  is  being  transmitted.  The  “minimum”  in  minimum-shift  keying 
refers  to  the  tone  spacing  t t/T,  which  is  smaller  than  that  used  in  most  FSK  systems. 


There  are  two  other  well-known  MSK-like  modulation  formats:  SFSK  (sinusoidal  FSK)  [2], 
described  by 


7 Tt  1  27 Tt 

2  T  4  Sm  T 


pit)  =  cos 

and  DSFSK  (doubly  sinusoidal  FSK)  [3],  described  by 
pit)  =  cos 


|/|  <  T, 


TTt_ 

IT 


1  .  2tt t  .  1  47 Tt 

—  sm  — —  +  — -  sm  — — 
3  T  24  T 


M  <  T. 


SFSK  and  DSFSK  are  not  discussed  further  except  for  comparison  of  spectral  rolloff  properties  with 
RSK  below. 


3.  THE  CLASS  OF  BASEBAND  PULSE  FUNCTIONS  OPTIMIZED 

Suppose  an  MSK-like  modulation  is  desired  using  the  implementation  strategy  outlined  in  Fig.  2. 
Further,  suppose  an  implementation  is  desired  with  the  Laplace  transform  of  fit),  denoted  by  Fis),  a 
rational  function  of  s.  (Throughout  this  report,  upper  and  lowercase  letters  refer  to  the  Laplace  and 
time  domains.)  This  would  be  necessary  if  the  arm  filters  were  to  be  implemented  using  conventional 
lumped-constant  filter  technology.  Are  there  any  obvious  requirements  on  Fis)  that  can  be  deduced 
from  the  constant-envelope  conditions  derived  above? 
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Equation  (11)  requires  that  pit )  be  supported  only  over  a  finite  interval.  This  implies  that  Pis) 
can  have  no  poles.  Since  Pis)  =  Gis)Fis),  Fis)  can  have  poles  only  where  Gis)  has  zeros  (located 
at  s  =  ±  j  kit  IT,  k  =  1, 2  =  . . .).  The  number  of  zeros  of  Fis)  can  be  determined  by  the  desired 
spectral-rolloff  properties  of  Pis).  Since  there  is  no  point  to  this  exercise  unless  the  spectrum  of  the 
resulting  baseband  pulse  function  rolls  off  faster  than  for  rectangular  QPSK,  Fis)  should  have  at  least 
one  more  pole  than  it  has  zeros.  Fis)  can  therefore  be  expressed  as 


Fis)  =  £ 

k>  1 


+  Pk 

s 2  +  ikir  /  T)2  ' 


where  the  cr’s  and  p’s  are  real  and  only  a  finite  number  of  terms  have  nonzero  cr  and  p.  Cascading  this 
Fis)  with  Gis)  effectively  substitutes  a  (smaller)  set  of  zeros  in  place  of  some  (finite)  set  of  the  zeros 
of  Gis).  The  most  obvious  choice  of  a  set  of  zeros  to  replace*  is  the  set  of  2 n  zeros  in  the  region  of 
the  s-plane  where  Gis)  is  generally  the  largest,  i.e.,  near  s  =  0,  effectively  making  the  above  sum 
range  from  k  —  1  to  k  —  n.  The  zeros  of  Fis)  can  be  further  restricted  by  requiring  <rk  =  0  in  the 
expression  above.  Necessary  for  computational  feasibility,  this  constrains  the  zeros  of  Fis)  to  occur  in 
pairs  on  the  axes  (real  or  imaginary)  symmetric  about  the  origin;  this  is  convenient  for  certain  popular 
realization  strategies.  The  following  results: 


Fis)  =  T  - — - 

£l  s2+  IknJW' 


(15) 


The  effect  on  pit)  of  this  last  restriction  is  to  force  it  to  be  symmetric  in  time,  i.e., 
pit0  +  t)  =  pit0-  t).  It  is  worth  noting  that  all  of  the  well-known  MSK-like  modulation  formats  are 
symmetric  in  time. 


The  class  of  functions  described  by  Eq.  (15)  can  be  described  more  conveniently  for  the  purposes 
of  this  report  as 


Fis)  = 


2t r2  "  ffjk 
P  £[  s2  +  ikir/T)2' 


(16) 


The  additional  constant  factor  and  the  change  from  pk  to  k2ak  will  allow  simpler  expressions  below. 

Although  the  use  of  arm  filters  described  by  Eq.  (16)  in  the  OQPSK  configuration  of  Fig.  2 
results  automatically  in  the  satisfaction  of  the  condition  expressed  in  Eq.  (11),  it  does  not  automatically 
result  in  the  satisfaction  of  Eq  (14).  An  optimization  criterion  will  therefore  be  chosen  in  the  following 
section  that  can  be  used  to  choose  the  coefficients  {fl^}  ih  Eq.  (16)  so  that  either  Eq.  (14)  Is  satisfied, 
or,  more  likely,  Eq.  (14)  is  nearly  satisfied.  There  are  two  other  constraints  that  need  to  be  satisfied. 
One  is  the  Nyquist  criterion  of  Eq.  (3),  necessary  to  prevent  intersymbol  interference.  The  remaining 
constraint,  necessary  to  ensure  that  the  OQPSK  signal  experiences  a  net  phase  change  of  90°  in  each  bit 
interval,  is  given  by  Eq.  (9).  Notice  that  any  function  pit)  satisfying  Eq.  (9)  satisfies  the  Nyquist  cri¬ 
terion  of  Eq.  (3)  automatically.  For  baseband  pulse  functions  satisfying  Eq.  (11)  the  additional  condi¬ 
tion 

pit0)  5*  0  (17) 

is  sufficient  (though  not  necessary)  to  ensure  that  Eq.  (9)  is  satisfied.  A  reasonable  strategy,  therefore, 
is  to  perform  the  optimization  first  and  then  check  to  see  whether  Eq.  (17)  (or,  if  not  Eq.  (17),  then 
Eq.  (9))  is  satisfied. 


*  Although  it  would  be  interesting  to  explore  optimum  arm-filter  functions  with  zeros  at  various  other  combinations  of  the  low- 
frequency  zeros  of  Gis),  such  a  task  is  beyond  the  scope  of  this  investigation, 
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It  is  convenient  to  express  the  optimization  criterion  in  terms  of  pit).  To  find  the  baseband  pulse 
function  pit)  implied  by  the  transfer  function  of  Eq.  (16),  the  impulse  response  corresponding  to  Eq. 
(16)  is  first  calculated  by  inverse  Laplace  transforming  Fis)  (in  a  region  of  convergence  Refs)  >  0): 

fit)  -  Uit)  £  kak  sin 

1  k- 1  1 

To  obtain  the  pulse  function  pit),  convolve  fit)  with  the  rectangular  driving  function  git).  Because 
fit)  is  periodic  with  period  IT  equal  to  the  width  of  the  gate  function  g  it),  the  convolution  will  equal 
zero  for  t  ^  2T: 

oo  t 

pit)  —  J  fir)  git  -  t)  dr  —  #(?)/  fir)  dr, 

—  oo  0 

or,  performing  the  integration, 

pit)  =  git)  £  2ak 

k- 1 

It  is  interesting  to  note  that  if  this  pulse  function  were  extended  by  letting  the  number  of  terms  n  in 
the  sum  approach  infinity,  it  would  become  a  Fourier  series  capable  of  representing  any  function  that 

(a)  was  real  with  finite-energy, 

(b)  had  support  only  over  the  interval  0  <  t  <  2  T,  and 

(c)  was  symmetric  with  respect  to  t  =  T,  i.e.,  piT  +  t)  =  piT  —  t). 


1  -  cos 


lent 


(18) 


4.  THE  OPTIMIZATION  CRITERION 

The  optimization  criterion  used  to  optimize  the  fa*.}  of  Eq.  (18)  must  be  one  which  tends  to 
represent  the  amplitude  fluctuations  in  the  OQPSK  signal  generated  with  the  baseband  pulse  function 
of  Eq.  (18).  To  be  able  to  deal  with  such  quantities  as  "the  average  fluctuation"  more  precisely,  it  is 
helpful  to  reformulate  the  complex  modulation  function  m  ii)  of  Eq.  (6)  as  a  random  process.  Let  the 
transmitted  data  symbols  /3m  be  independent  random  variables  with  the  two  possible  values  equiprob- 
able,  and  interpret  the  summation  over  the  index  m  as  ranging  from  m  =  -  »  to  m  =  ~.  Make  this 
random  process  stationary  by  inserting  a  random  delay  u,  uniformly  distributed  from  0  to  T.  These 
changes  produce 


mR  it)  =  ~  mT  —  «);  W~U(0,77). 

m 


(19) 


The  subscript  on  the  function  mR  ( t )  indicates  a  random  process. 


Given  this  formulation,  one  obvious  strategy  is  to  choose  the  { ak }  of  Eq.  (18)  to  minimize  the 
variance  of  \mR(t)\  while  holding  the  mean  of  |m*G)|  constant.  Unfortunately,  the  mean  of  \mR(t)\ 
is  difficult  to  calculate,  A  slightly  modified  strategy  proves  easier:  choose  the  { ak }  of  Eq.  (18)  to 
minimize  the  variance  of  \mR(t)\2  while  holding  the  mean  of  \mR(t)\2  constant.  The  extra  squaring 
has  the  added  benefit  of  increasing  the  "penalty”  associated  with  large  envelope  fluctuations.  Appendix 
A  shows  that  the  mean  of  \mR  (t)  I2  is 


n 


X  °k 

k- 1 


+  4  £  al. 


(20) 
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and  the  variance  of  \mR  it)  |2  is 


n  n  n 


Var{|w*0)|2}  =  £  £  £  i,akaiam ar 

/c-1  /« 1  m-1  r—  1 


27  even  (A:)  even  (/)  8^_  ,■ 

-  24  8*_,8m_r 

-  2 6 even (k)  i8k+Hm+8k-l+m+8k-Hm) 
2  even  (&+/)  (  8<r+/+m— »«+/• 

&k+l—  m—  r^~  &k—  i+m+r 
~^~^k—  i+m-  ®lt-  i—  m+r 

+  8* /-m-f) 


(21) 


where  the  function  even  ( )  is  defined  to  be  unity  if  its  integer  argument  is  even  and  zero  otherwise. 
Equations  (20)  and  (21)  will  be  evaluated  separately  below  for  the  4th-,  6th-,  and  8th-order  arm  filters 
discussed  in  the  next  section. 


Optimization  Constraints  to  Control  Spectral  Rolloff 

Since  Fis)  in  Eq.  (16)  has,  in  general,  two  more  poles  than  zeros,  its  magnitude  rolls  off  asymp¬ 
totically  with  frequency  at  12  dB/octave.  Since  Gis )  rolls  off  at  6  dB/octave,  Pis)  rolls  off  at  18 
dB/octave,  faster  than  the  12  dB/octave  of  MSK  but  slower  than  the  24  dB/octave  of  SFSK.  If  faster 
spectral  rolloff  than  18  dB/octave  is  needed,  the  {ak}  will  have  to  be  further  constrained  during  optimi¬ 
zation  to  force  F(s )  to  have  fewer  zeros. 


To  see  what  these  constraints  are,  it  is  helpful  to  express  the  arm-filter  transfer  function  Fis)  of 
Eq.  (16)  as  a  ratio.  Giving  F(s)  a  common  denominator  produces 


Fis )  = 


27T2 

T 2 


t  k2ak  fls2+ihr/T)2 

k-l  /-I 

i*k 

s2  +  ( kn/T )2 

k-l 


(22) 


Before  trying  to  put  the  numerator  of  this  expression  in  the  form  of  a  single  polynomial  in  s,  consider 
expanding  the  simpler  product, 


IT(s  +  • 

k-l 


This  expands  into  a  nth-degree  polynomial  in  s,  where  the  coefficient  of  s"  m  is 


b, 


*1**2' 


■A  •  •  A 


with  defined  to  be  the  set  of  all  combinations  of  m  of  the  integers  between  1  and  n.  When  m  is 
zero,  'Fm  is  defined  to  contain  only  a  single,  empty  combination.  This  makes  the  product  of  the  {b/jt} 
empty,  and  thus  equal  to  the  identity  element  for  multiplication:  unity.  Using  this  approach,  the  coef¬ 
ficient  of  s2(n~  m)  in  the  expansion  of 

]Js2  +  Utt/T)2 

k-l 
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becomes 


/2_2 
l\  7 T 


('1.'2 . 


;  2  2 

1 2  7T 


jQ.  y2 


j  2  __  2 

T 2 


or 


7T 


2m 


j"2m 


('1'2  •••  4,): 


. 


To  modify  this  to  take  account  of  the  /  ^  k  condition  in  the  product  in  the  numerator  of  Eq.  (22), 
replace  the  set  of  combinations  T'"  with  a  modified  set  0">*  defined  to  be  the  set  of  all  combinations 
of  m  of  the  integers  between  1  and  n  excluding  k.  The  expression  above  will  then  give  the  coefficient 
of  s2(»-m-D  in  the  expansion  of  the  product.  Again,  ©o’’*  is  the  set  containing  only  the  empty  combi¬ 
nation.  Equation  (22)  can  now  be  rewritten  as 


F(s)  = 


2tt 2 
T2 


n- 1 

I5 

m-0 


2(n-  m—  \)JL 


2m 


j^2n 


X  k2ak 

k-\ 


X  (hh  ■■■  4.) 2 


ns2+  (kir /T)2 

k- 1 


(23) 


To  obtain  faster  spectral  rolloff  of  the  baseband  pulse  function  than  18  dB/ octave  requires  con¬ 
straining  the  { ak }  in  the  optimization  process  so  that  the  degree  of  the  numerator  of  the  arm-filter 
transfer  function  of  Eq.  (23)  is  reduced  to  less  than  2 n  -  2.  To  reduce  the  degree  of  the  numerator  by 
two  to  gain  an  additional  12  dB/ octave  requires  setting  the  coefficient  of  the  highest  degree  term  in  the 
numerator  to  zero, 

£*2afc-  0.  (24) 

A:— 1 

Notice  that  the  rightmost  sum  in  Eq.  (23)  contained  only  an  empty  product.  To  reduce  the  degree  of 
the  numerator  of  Eq.  (23)  by  four  to  gain  24  dB/octave  requires  Eq.  (24)  and  setting  the  next  coeffi¬ 
cient  to  zero  by  requiring 

0.  (25) 

Jt-1  1-1 

i^k 

The  only  limit  to  the  number  of  such  constraints  that  can  be  applied  is  that  at  least  one  degree  of 
freedom  among  the  {ak}  must  remain  for  optimization. 

5.  THE  RSK  MODULATION  FORMATS 

The  optimization  of  the  { ak }  of  Eq.  (16)  (and  Eq.  (18))  according  to  the  optimization  criterion 
developed  above  is  discussed  separately  below  for  4th-,  6th-,  and  8th-order  arm  filters.  In  each  case, 
the  following  information  is  given: 

•  A  name  for  each  format,  RSK(m.n),  where  m  and  n  indicate  the  numbers  of  zero  pairs  and  pole 
pairs  in  the  rational  transfer  function. 

•  The  optimum  coefficients  {ak}  for  use  in  Eqs.  (16)  and  (18). 

•  A  plot  of  the  baseband  pulse  function  pit)  given  in  Eq.  (18). 
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•  A  plot  of  one  period  of  the  magnitude  of  the  complex-baseband  OQPSK  signal  of  Eq.  (6).  This  is 
given  by  Eq.  (13)  with  t0  =  T. 


The  normalized  zero  locations  { zk }  and  gain  y  for  the  following  factored  equivalent  of  Eq.  (16): 


m 

n 

Ar— 1 

li 

zk 

2 

+ 1 

n 

n 

k- 1 

Ts 

kir 

1 

+ 1 

F(s )  =  y  - 7 - -  .  (26) 


A  plot  of  the  spectral  density  of  the  OQPSK  process  of  Eq.  (19),  shown  in  Appendix  B  to  be 


n 


S  •(ftj’)  —  4y2 

mRmR 


sin'2  ft  T  k- 


ft? 


n 

n 

k~  1 


1  - 


ft 


*k 


1  - 


ft 


k  v 


(27) 


The  spectral  densities  are  shown  in  a  form  in  which  the  spectral  level  is  given  as  a  function  of  a 
dimensionless  frequency  variable  ft  T  =  ft  T.  This  has  the  effect  of  making  the  density  dimension¬ 
less  as  well.* 


All  of  the  above  items  except  the  plots  are  repeated  in  tabular  form  in  the  Summary. 

Fourth-Order  Arm  Filters 


Fourth-order  arm  Filters  are  obtained  by  setting  n  =  2  in  Eqs,  (16)  and  (18).  Evaluating  Eq,  (20) 
with  n  =  2  yields  the  mean  squared  amplitude 

E  (I/MrO)!2}  =  12a2  +  160^3  +  12a2.  (28) 

Similarly,  evaluating  Eq.  (21)  with  n  =  2  yields  the  variance  of  the  squared  amplitude 

Var{|aiR  (r) |2}  =  136a4  +  256aia|  +  64 a?  a2  —  64a3  a2  +  8a4.  (29) 


The  optimization  problem  can  be  formulated  using  a  Lagrange  multiplier.  Define  the  Lagrangian 

L  =  Var{K(r)|2}  -/x(£{|wR(r)|2}  -  1)  (30) 

or,  substituting  from  Eqs.  (29)  and  (28), 

L  =  136a4  +  256aia?  +  64a2  a22  —  64a3  a2  +  8a4 
+  fi  ( 12  a22  +  16aia2  +  12a2  —  1). 

The  variable  fj.  is  a  Lagrange  multiplier  introduced  to  force  the  factor  that  it  multiplies  to  equal  zero  in 
the  optimization  solution.  The  optimization  now  proceeds  by  setting  the  gradient  of  the  Lagrangian 
(with  respect  to  the  three  variables  ah  a2,  and  /x)  to  zero  and  solving  for  the  stationary  points  of  the 
Lagrangian.  Computing  the  gradient  and  equating  to  zero  produces 

0  =  256a23  +  128a,a22  -  192a,2  a2  +  32«i3  +  fx(\6a2  +  24a,) 

0=  544a3  +  768a,  a22  +  128a?  a2  -  64a}3  +/*(24a2+  16a!) 

0  =  12a?  +  16aia2  +  12a2  —  1 . 

*Of  course,  if  the  pulse  function  were  to  be  given  units,  the  corresponding  normalized  spectral  densities  would  have  units  as  well. 
For  example,  if  p(t)  were  expressed  in  volts,  the  corresponding  normalized  spectral  density  would  have  units  of  volts2, 
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Solving  this  set  of  simultaneous  equations  numerically  produced  the  stationary  points  given  in  Table  1. 
The  negative  of  each  solution  shown  is  also  a  stationary  point  (of  the  same  type  and  with  the  same  vari¬ 
ance)  due  to  the  symmetry  of  the  problem.  As  the  choice  of  the  sign  of  p  ( t )  is  irrelevant  in  most  prac¬ 
tical  applications,  only  one  of  each  such  pair  of  solutions  is  given. 

Table  1  —  Lagrangian  Stationary  Points 
with  n  —  2 


a\ 

a2 

Var{  mR(t) |2} 

0.25143 

0.31903 

0.38542 

0.05863 

0.05195 

-0.37636 

-0.22858 

-0.32443 

5.85147  x  10~5 
0.16243 

0.70375 

1.02166 

A  discriminant  test  on  the  constraint  equation 

12  a2  +  16  a  i  #2  +  1 2  a  i2  —  1 

shows  that  it  describes  an  ellipse  in  the  (aua2)  plane.  Since  it  describes  a  closed  line,  the  global 
minimum  of  the  variance  on  that  line  is  just  the  stationary  point  for  which  the  variance  is  smallest,  i.e., 

ai  =  0.25143,  a2  =  0.05195. 

The  baseband  pulse  function  obtained  by  substituting  these  values  into  Eq.  (18)  with  n  «  2  is  plotted 
in  Fig.  3.  Since  it  is  nonzero  in  the  center  as  required  by  Eq.  (17),  it  satisfies  the  90°  net-phase-change 
criterion  of  Eq.  (9)  and  the  Nyquist  (zero-intersymbol-interference)  criterion  of  Eq.  (3)  automatically. 
Figure  4  shows  the  i(t)  and  q  it)  waveforms  of  Fig.  2  that  would  result  from  using  the  arm  filters  that 
correspond  to  this  baseband  pulse  function. 


Fig.  3  —  Baseband  pulse  function  from  optimum  fourth-order  arm  Filter 


Figure  5  shows  one  period  of  the  magnitude  (calculated  using  Eq.  (13))  cf  the  complex-baseband 
OQPSK  signal  of  Eq.  (6)  using  the  baseband  pulse  function  of  Fig.  3.  This  shows  that  the  square  root 
of  the  sum  of  the  squares  of  the  i  and  q  waveforms  of  Fig.  4  would  fluctuate  less  than  1%  from  its 
mean.  The  amplitude  variations  are  so  small  as  to  be  practically  negligible. 

If  the  transfer  function  of  the  arm  filter  is  put  into  the  factored  form  of  Eq.  (26),  there  is  one  pair 
of  imaginary  zeros  (m  —  1)  and  the  gain  constant  and  zero  location  are  given  by 

y  —  0.606768,  z\  —  1.62557  tt. 
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Fig.  4  —  Typical  i  and  q  waveforms  using  the  fourth-order  arm  filters 


Fig.  5  —  One  period  of  the  envelope  of  the  fourth-order  RSK  signal 


The  spectral  density  of  the  fourth-order  RSK  process  is  shown  in  Fig.  6  along  with  the  spectral 
density  of  the  MSK  and  SFSK  processes,  for  comparison.  All  three  densities  are  normalized  to  unit 
total  power.  The  RSK  density  appears  to  be  between  the  MSK  and  SFSK  densities  in  several  ways: 
spectral  density  at  zero  frequency,  width  of  main  spectral  lobe,  height  of  First  sidelobe,  and  asymptotic 
rolloff  rate.  Thus,  the  fourth-order  RSK  modulation  format  can  be  used  where  a  compromise  between 
the  spectral  properties  of  MSK  and  SFSK  is  needed. 


2w 

Fig.  6  —  Spectral  densities  of  fourth-order  RSK  (solid), 
MSK  (dashed),  and  SFSK  (chain-dashed) 
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Sixth-Order  Arm  Filters 

Sixth-order  arm  filters  are  obtained  by  setting  n  =  3  in  Eqs.  (16)  and  (18).  Evaluating  Eq.  (20) 
with  n  =  3  yields  the  mean  squared  amplitude 

E{\mR  it)  |2)  =  12a32  +  4a3  (4a3  +  4a2)  +  16a2a3  +  12  a2  4-  12a2. 

Similarly,  evaluating  Eq.  (21)  with  n  =  3  yields  the  variance  of  the  squared  amplitude 

Var{|oi*  ( t ) |2)  =  8a3  +  128a2  a3  —  128a!a2a32  -I-  64a2  a3  +  256a2  a3  +  160a!a2  a3 

—  192a2  a2a3  +  32a2  a3  +  136a2  +  256a!a2  +  64a2  a2  —  64a2  a2  +  8a4. 


The  optimization  of  this  variance  with  only  a  mean-squared  amplitude  constraint  would  lead  to  a 
pulse  pit )  with  an  asymptotic  spectral  rolloff  rate  no  better  than  that  associated  with  the  fourth-order 
problem  solved  earlier.  The  extra  degree  of  freedom  inherent  in  the  sixth-order  problem  would  be 
expected  to  lead  to  a  solution  with  lower  variance,  but,  since  the  variance  obtained  by  the  simpler 
fourth-order  solution  is  adequate  for  any  practical  application,  there  is  little  point  in  looking  for  a  higher 
order  filter  function  with  the  same  rolloff  rate.  To  obtain  a  higher  rolloff  rate,  the  constraint  given  by 
Eq.  (24)  can  be  used  to  force  the  sixth-order  transfer  function  F(s)  to  have  four  more  poles  than 
zeros.  Solving  Eq.  (24)  (with  n  =  3)  for  al  gives 

a\  —  —  9  a3  —  4a2. 

Substitution  of  this  expression  for  ax  into  the  expressions  just  given  for  the  mean  and  variance  of  the 
squared  amplitude  gives  the  constrained  mean 

E  it)  I2)  =  840 a3  +  672a2a3  +  140 a2 , 

and  the  constrained  variance 

Var{|m/{  (r)|2}  =  34352 a3  +  99072a2a33  +  102176a22  a32  +  42880a23  a3  +  6280a24. 


The  optimization  problem  can  now  be  formulated  using  a  Lagrange  multiplier.  The  use  of  the 
Lagrangian  form  given  in  Eq.  (30)  with  the  constrained  mean  and  variance  above  gives 

L  =  34352a3  +  99072a2a33  +  102176a22  a32  +  42880a23  a3  +  6280a24 
T  p  i 840 a 2  672a2a3  -I-  140 a 2  —  1 ) . 

Computing  the  gradient  of  this  Lagrangian  and  equating  it  to  zero  produces  three  simultaneous  non¬ 
linear  equations  in  three  unknowns: 

0  =  99072a3  +  204352a2a32  +  128640a22  a3  +  25120a23 
+  /z  (672 a3  +  280a2) 

0  =  137408 a3  +  297216a2a32  +  204352a22  a3  +  42880a23 
+  p,  (1680a3  +  672a2) 

0  =  840a3  +  672a2a3  +  140a2  —  1. 

Solving  this  set  of  simultaneous  equations  numerically  produced  the  stationary  points  given  in  Table  2. 
Again,  due  to  the  symmetry  inherent  in  the  problem,  the  negative  of  each  solution  shown  is  also  a  sta¬ 
tionary  point. 
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Table  2  —  Lagrangian  Stationary  Points  with 
_ n  =  3  and  a1  Constrained 


a\ 

a2 

«3 

Vardm^  (?)  I2} 

0.29031 

-0.18456 

0.04338 

0.34040 

0.04595 

0.42203 

0.34356 

-0.29482 

-0.05268 

-0.16706 

-0.15751 

0.09321 

6.56782  x  10-3 
0.47046 

0.94788 

1.12266 

A  discriminant  test  on  the  constraint  equation 

840a32  +  672a2a3  +  140a2  =  1 

shows  that  it  describes  an  ellipse  in  the  ( a2,a3 )  plane.  Since  it  describes  a  closed  line,  the  global 
minimum  of  the  variance  on  that  line  is  just  the  stationary  point  for  which  the  variance  is  smallest,  i.e., 

ax  =  0.29031,  a2  =  0.04595,  a3  =  -0.05268. 

The  baseband  pulse  function  obtained  by  substituting  these  values  into  Eq.  (18)  with  n  =  3  is  plotted 
in  Fig.  7.  Since  it  is  nonzero  in  the  center  as  required  by  Eq.  (17),  it  satisfies  the  90°  net-phase-change 
criterion  of  Eq.  (9)  and  the  Nyquist  (zero-intersymbol-interference)  criterion  of  Eq.  (3)  automatically. 


The  pronounced  dip  in  the  middle  of  the  pulse  function  shown  in  Fig.  7  is  consistent  with  the 
relativley  large  (squared-amplitude)  variance  implied  by  Table  2.  This  effect  is  shown  more  clearly  in 
Fig.  8,  which  shows  one  period  of  the  magnitude  (calculated  using  Eq.  (13))  of  the  complex-baseband 
OQPSK  signal  of  Eq.  (6)  using  the  baseband  pulse  function  of  Fig.  7.  While  this  high  envelope  ripple 
may  be  intolerable  in  some  applications,  it  may  be  possible  to  tolerate  it  in  others  in  order  to  obtain  the 
associated  spectral  characteristics  given  in  Fig.  9. 


The  spectral  density  of  the  sixth-order  RSK  process  in  Fig.  9  is  shown  along  with  the  spectral  den¬ 
sity  of  the  SFSK  and  DSFSK  processes,  for  comparison.  All  three  densities  are  normalized  to  unit  total 
power.  The  RSK  density  appears  to  be  between  the  SFSK  and  DSFSK  densities  in  several  ways:  spec¬ 
tral  density  at  zero  frequency,  width  of  main  spectral  lobe,  height  of  first  sidelobe,  and  asymptotic  roll¬ 
off  rate.  The  RSK  spectrum  is,  however,  clearly  superior  to  the  DSFSK  spectrum  for  the  first  few 
sidelobes.  Thus,  while  the  DSFSK  spectrum  is  generally  thought  so  poor  as  to  preclude  the  implemen¬ 
tation  of  DSFSK,  the  RSK(2,6)  format  may  be  useful  where  a  higher  rolloff  rate  is  needed  than  is  pro¬ 
vided  by  SFSK  and  where  the  substantial  ripple  of  RSK(2,6)  is  tolerable. 


If  the  transfer  function  of  the  arm  filter  is  put  into  the  form  given  by  Eq.  (26),  there  is  one  pair 
of  imaginary  zeros  {m  =  1)  and  the  gain  constant  and  zero  location  are  given  by 

y  =  0.56716,  z,  =  1.77465 7r. 
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Fig.  8  —  One  period  of  the  envelope  of  the  sixth-order  RSK  signal 


2r 

Fig.  9  —  Spectral  densities  of  sixth-order  RSK  (solid), 
SFSK  (dashed),  and  DSFSK  (chain-dashed) 


Eighth-Order  Arm  Filters 

To  obtain  the  30  dB/octave  asymptotic  spectral  rolloff  of  the  previous  section  without  the  associ¬ 
ated  envelope  ripple  requires  adding  another  degree  of  freedom  to  the  optimization  by  moving  from  a 
sixth-order  arm  filter  to  an  eighth-order  arm  filter.  Eighth-order  arm  filters  are  obtained  by  setting 
n  =  4  in  Eqs.  (16)  and  (18).  Evaluating  Eq.  (20)  with  n  =  4  yields  the  mean  squared  amplitude 

E  { |  mj?  (t)  |2)  =  12a4  +  16a3a4  +  16a2a4  +  16aia4  +  12 a 2  +  16a2a3  +  16# ^3  +  12  a2 

+  16c1o2  +  12a2. 

Similarly,  evaluating  Eq.  (21)  with  n  —  4  yields  the  variance  of  the  squared  amplitude 

Var{|w/{ (r) |2}  =  136a4  +  256a3a4  +  256 a2a4  +  25 6aia4  +  128 <23  a\  +  25 6a2a3a4 
+  128a!a3a4  +  128a2  a4  +  256aja2a4  +  128a2  a4  +  32 a2a3  a4 
—  128aja32  a4  +  64a2  a3a4  —  192a!a2a3a4  —  128a2  a3a4  +  64a2  a4 
+  64a]a2  a4  -  32a2  a2a4  +  8a3  +  128a2  a3  —  128a!a2a3  +  64a2  a3 
+  256 a2  a3  +  160aia2  a3  —  192a2  a2a3  +  32a2  a3  4-  136a2  4-  256aia| 

4-  64a2  a 2  —  64a2  a2  +  8af. 


15 


J.  O.  COLEMAN 


The  optimization  of  this  variance  with  only  a  mean-squared  amplitude  constraint  would  lead  to  a 
pulse  pit)  with  an  asymptotic  spectral  rolloff  rate  no  better  than  that  associated  with  the  fourth-order 
problem  solved  earlier.  To  obtain  a  higher  rolloff  rate,  the  constraint  given  by  Eq.  (24)  can  be  used  to 
force  the  eighth-order  transfer  function  Fis)  to  have  four  more  poles  than  zeros.  Solving  Eq.  (24) 
(with  «  =  4)  for  a\  gives 


a  i=  —  16a4  —  9a3  —  4a2. 

Substitution  of  this  expression  for  ax  into  the  expressions  just  given  for  the  mean  and  variance  of  the 
squared  amplitude  gives  the  constrained  mean 

E  {| mu  (r)|2)  =  2828 <74  +  3072 a3a4  +  1232a3a4  +  840 a2  +  £ >72 <22 a 3  +  140a2 
and  the  constrained  variance 

Var{|m*(f)|2}  -  553096a44  +  1048576a3a43  +  789760a2a43  +  765056a32a42 

+  1163776a2a3a4  +  405632a22  a4  +  258048a3  a4  +  581312a2a3  a4 
-I-  404992a22  a3a4  +  85312a23  a4  +  34352a4  +  99072 a2a3  +  102176a2  a3 
+  42880a23a3  +  6280a24. 


In  theory,  the  optimization  problem  can  now  be  formulated  using  a  Lagrange  multiplier  as  was 
done  earlier  when  solving  for  the  fourth-  and  sixth-order  arm  filters.  In  practice,  however,  the  result¬ 
ing  set  of  four  simultaneous  nonlinear  equations  in  four  unknowns  proved  too  difficult  to  solve  using 
the  computational  facilities  at  hand. 

As  an  alternative,  an  iterative  search  program  was  used  to  search  the  (a2,  a3)  space  for  a 
minimum  of  the  constrained  variance  expression  above.  For  each  (a2,  a3)  point,  a4  was  calculated  by 
evaluating  the  constrained  mean  (given  above)  at  the  (a2,  a3)  point  in  question,  setting  the  resulting 
expression  to  unity,  and  solving  the  resulting  quadratic  equation  in  a4.  Each  point  in  (a2,  a3)  space 
thus  provided  two  points  in  (a3,  a3,  a4)  space  at  which  to  evaluate  the  constrained  variance. 

From  a  starting  point  at  which  a3  and  a3  were  equal  to  the  values  given  earlier  as  the  optimums 
for  the  sixth-order  arm  filter,  the  search  rapidly  converged  on  a  minimum,  shown  in  Table  3.  The 
corresponding  variance  was  1.25202  x  10~6,  smaller  even  than  the  fourth  order  optimum.  This  is  not 
surprising,  in  view  of  the  extra  degree  of  freedom  here. 


Table  3  —  Minimum-Variance  Point 
with  n  =  4  and  a\  Constrained 


k 

1 

0.28130 

2 

0.05200 

3 

-0.03141 

4 

-0.01291 

Due  to  the  nature  of  the  method  employed  to  find  this  minimum,  there  is  no  way  to  be  certain  it 
is  a  global  minimum.  This  is  not  important,  however,  because  this  solution  does  have  the  properties 
desired.  The  baseband  pulse  function  pit)  obtained  by  substituting  the  values  of  Table  3  into  Eq.  (18) 
with  n  =  4  is  plotted  in  Fig.  10.  Once  again,  since  it  is  nonzero  in  the  center  as  required  by  Eq.  (17),  it 
satisfies  the  90°  net-phase-change  criterion  of  Eq.  (9)  and  the  Nyquist  (zero-intersymboMnterference) 
criterion  of  Eq.  (3)  automatically. 


16 


NRL  REPORT  8871 


Fig.  10  —  Baseband  pulse  function  from  optimum  eighth-order  arm  filter 

Figure  11  shows  one  period  of  the  magnitude  (calculated  using  Eq.  (13))  of  the  complex- 
baseband  OQPSK  signal  of  Eq.  (6)  using  the  baseband  pulse  function  of  Fig.  10.  The  amplitude  varia¬ 
tions  are  even  smaller  than  those  shown  in  Fig.  5  for  the  fourth-order  case. 

If  the  transfer  function  of  the  arm  filter  is  put  into  the  form  given  by  Eq.  (26),  there  are  two  pairs 
of  imaginary  zeros  ( m  =  2)  and  the  gain  constant  and  zero  locations  are  given  by 

y  =  0.57795,  z,  -  1.72733  tt,  and  z2  =  3.43183  n. 

I  .01  I —  - *1 


1 


Fig.  11  —  One  period  of  the  envelope  of  the  eighth-order  RSK  signal 


The  spectral  density  of  the  eighth-order  RSK  process  is  shown  in  Fig.  12  along  with  the  spectral 
density  of  the  SFSK  and  DSFSK  processes,  for  comparison.  All  three  densities  are  normalized  to  unit 
total  power.  Like  the  sixth-order  RSK  spectrum  shown  in  Fig.  9,  the  eighth-order  RSK  spectrum 
appears  to  be  between  the  SFSK  and  DSFSK  densities  in  several  ways:  spectral  density  at  zero  fre¬ 
quency,  width  of  main  spectral  lobe,  height  of  first  sidelobe,  and  asymptotic  rolloff  rate.  The  RSK 
spectrum  is,  however,  clearly  superior  to  the  DSFSK  spectrum  for  the  first  few  sidelobes.  Figure  13 
illustrates  the  tradeoff  between  the  sixth-  and  eighth-order  RSK  formats  by  showing  their  spectra 
together  for  comparison.  Clearly,  the  more  constant  envelope  attained  in  the  eighth-order  case  was  at 
the  expense  of  desired  spectral  qualities.  The  eighth-order  RSK  spectrum  has  a  substantially  higher 
second  sidelobe,  as  well  as  somewhat  higher  sidelobes  in  the  high-frequency  rolloff  region.  The 
eighth-order  RSK  format  is  therefore  preferred  over  the  sixth-order  format  only  when  the  much  lower 
envelope  ripple  justifies  both  the  added  complexity  and  a  somewhat  poorer  spectrum. 

6.  SUMMARY 

Three  new  baseband  pulse  functions  were  derived  for  OQPSK  signalling.  I  describe  them  as 
“almost  MSK-like”  because  the  OQPSK  signals  generated  using  these  pulse  functions  (a)  change  phase 
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AT 


Fig.  12  -  Spectral  densities  of  eighth-order  RSK  (solid), 
SFSK  (dashed),  and  DSFSK  (chain-dashed) 


Fig.  13  —  Spectral  densities  of  sixth-order  RSK  (solid) 
and  eighth-order  RSK  (dashed) 


by  exactly  90°  in  a  bit  interval  and  (b)  have  almost  constant  amplitude.  In  addition,  the  new  formats 
have  asymptotic  spectral  rolloff  rates  greater  than  the  12  dB/octave  of  MSK.  These  properties  in 
OQPSK  modulation  types  generally  result  in  robust  waveforms  better  able  to  tolerate  the  nonlinearitie* 
in  high-power  traveling-wave  tube  (TWT)  amplifiers  than  conventional  modulation  types.  Thus,  these 
new  modulation  formats  should  prove  especially  useful  where  TWT  operation  at  saturation  is  desired 
(and  the  wider  central  spectral  lobe  can  be  tolerated),  such  as  in  pulse-compression  radar  systems  and 
in  certain  satellite  data-communication  systems. 
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The  mathematical  form  of  the  new  pulse  functions  is  a  sum  of  raised  cosines: 

,  0<  t < 2T ; 


pit)  - 

2  ak 

k- 1 

1 

O 

o 

C/5 

0  elsewhere. 

The  form  itself  suggests  implementation  of  these  formats  by  physical  generation  of  the  raised-cosine 
waveforms  for  summing  into  two  trains  (in-phase  and  quadrature  baseband  components)  of  shaped 
pulses  that  can  be  polarity-switched  by  the  data  streams.  This  is  similar  to  a  commonly  used  strategy 
for  generating  MSK. 


Another  method  of  generating  these  modulation  formats  is  suggested  (from  Eqs. 
(26)  and  Fig.  2)  by  writing  the  Laplace  transform  of  the  baseband  pulse  function  p  (/)  in 
form: 


F{s)  = 


(1  -  e~s2T ) 


m 

y  n 

£-1 

Ts_ 

zk 

2 

+  1 

n 

n 

Ts 

ICTT 

1 

+  1 

(5),  (8),  and 
the  following 


This  is  the  Laplace  transform  of  the  rectangular  pulse  function  (of  rectangular  OQPSK)  multiplied  by  a 
rational  transfer  function.  This  transfer  function  can  be  approximately  implemented  (putting  the  poles 
as  close  to  the  imaginary  axis  as  possible)  in  simple  arm  filters  and  used  to  shape  NRZ  data  waveforms 
in  the  /  and  q  channels.  I  have  called  these  new  modulation  formats  rational-shift  keying,  emphasizing 
both  the  similarity  to  MSK  and  the  rational  transfer  functions  in  the  Laplace  transform  of  the  pulse  for¬ 
mats.  Each  format  is  denoted  by  RSK(m,n),  where  m  and  n  indicate  the  numbers  of  zero  pairs  and 
pole  pairs  in  the  rational  transfer  function. 


The  parameters  n,  ak,y,  m,  and  zk  in  the  above  equations  are  shown  in  Table  4  and  were  chosen 
(for  each  combination  of  m  and  n)  to  minimize  the  variance  of  the  amplitude  squared  of  a  unit-power 
baseband  OQPSK  random  process.  Along  with  the  RSK  parameters,  Table  4  gives  this  associated  vari¬ 
ance,  the  peak-to-peak  amplitude  ripple,  the  amount  the  first  spectral  sidelobe  is  down  from  the  peak, 
and  the  asymptotic  spectral  rolloff  rate. 


Table  4  —  Summary  of  RSK  Parameters  and  Characteristics 


Parameter 

RSK  (2,4) 

RSK  (2,6) 

RSK(4,8) 

a\ 

02 

04 

0.25143 

0.05195 

0.29031 

0.04595 

—0.05268 

0.28130 

0.05200 

-0.03141 

-0.01291 

y 

Z\ 

*2 

0.606768 

1.62557tt 

0.56716 

1.77465tt 

0.57795 

1.72733tt 

3.43183tt 

Vardwi/jWl2} 

Peak-to-peak  envelope  fluctuation 

5.85147  x  10~5 
0.0111 

6.56782  x  10“3 
0.1237 

1.25202  x  10~6 

0.0019 

Spectral  rolloff  rate 
Peak-to-first-sidelobe  ratio 

18  dB/octave 

21.0  dB 

30  dB/octave 

17.12  dB 

30  dB/octave 

18.04  dB 

OQPSK  spectral  density 

Baseband  pulse  function 

Envelope 

Figure  6 

Figure  3 

Figure  5 

Figure  9 

Figure  7 

Figure  8 

Figure  12 

Figure  10 

Figure  11 

19 


J.  O.  COLEMAN 


7.  REFERENCES 

1.  S.  Pasupathy,  "Minimum  Shift  Keying:  A  Spectrally  Efficient  Modulation,"  IEEE  Commurt.  Mag. 
17(4),  14-22  (July  1979). 

2.  F.  Amoroso,  "Pulse  and  Spectrum  Manipulation  in  the  Minimum  (Frequency)  Shift  Keying 
(MSK)  Format,"  IEEE  Trans.  Commun.  COM-24(3),  381-384  (March  1976). 

3.  B.  Bazin,  "A  Class  of  MSK  Baseband  Pulse  Formats  with  Sharp  Spectral  Roll-Off,"  IEEE  Trans. 
Commun.  COM-27  (5),  826-829  (May  1979). 


20 


Appendix  A 

THE  MEAN  AND  VARIANCE  OF  THE  SQUARED  AMPLITUDE 

In  this  appendix,  I  calculate  the  mean  and  variance  of  the  squared  amplitude  of  the  OQPSK  pro¬ 
cess  described  by  Eqs.  (18)  and  (19).  I  calculate  the  variance  by  calculating  the  second  and  fourth 
moments  of  the  amplitude  of  the  process  and  then  using 

Var{|m*U)l2}  =  E  dm*  it)  |4}  -  {E{\mR(t)\2))2,  (Al) 

to  obtain  the  variance  from  those  second  and  fourth  moments. 


The  Fourth  Moment  of  the  Amplitude 

I  first  calculate  the  fourth  moment  of  the  amplitude  of  the  OQPSK  process  described  by  Eq.  (19) 
with  no  restrictions*  on  the  real  baseband  pulse  function  pit).  The  pit)  of  Eq.  (18)  is  considered  later 
as  a  special  case  of  the  general  result. 

The  General  Result:  Real  Baseband  Pulse  Functions 


The  integration  over  the  uniform  random  variable  u  in  Eq.  (19)  is  made  explicit  when  the  desired 
moment  is  formulated  in  terms  of  a  conditional  expectation: 

T 


E  {\mR  (f)  |4}  =  ~  j  £{|/M/{0)l4 


u)  du 


(A2) 


Now,  calculate  that  conditional  expectation  by  substituting  for  mR{t)  from  Eq.  (19)  and  factoring  all 
the  nonrandom  variables  (including  «,  in  this  context)  out  of  the  expectation: 


E{\mR  (r)  |4 


wi¬ 


ll  I  it  £{/3,/MWU//rmrfl 

i  k  m  n 


•  pit  —  iT  —  u)  pit  —  kT  —  u)  pit  —  mT  —  u)  pit  —  nT  —  u)  ].  (A3) 

The  baseband  pulse  function  pit)  has  been  assumed  real.  Consider  the  remaining  expectation 
E  {/3//3/f/3m/3n}.  Because  the  expectation  of  a  product  of  independent  random  variables  is  equal  to  the 
product  of  the  expectations  of  those  variables,  and  because  the  expectations  of  each  of  the  {fik}  is  zero, 
E  {fijfikfimfin}  will  be  zero  whenever  at  least  one  of  the  indices  /,  k ,  m ,  and  n  is  different  from  all  of 
the  others.  In  the  complementary  case  in  which  each  index  is  equal  to  at  least  one  of  the  others,  each 
ft  appearing  in  the  product  appears  twice,  removing  the  sign  information.  Therefore,  using  symbols  A 
for  logical  conjunction  (“and”)  and  V  for  logical  disjunction  (“or”), 

V  /  = k A  m  =  n  , 

V  /  =  m  A  k—  n  , 

V  /  =  n  A  k—  m\ 
otherwise  . 


EifiiPkPmPn)  = 


*Other  than  those  necessary  to  assure  the  convergence  of  various  sums  and  integrals;  convergence  is  not  discussed  explicitly  but 
is  assured  for  pit)  of  practical  interest. 
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Noting  that  the  three  cases  where  the  expectation  equals  unity  overlap  only  where  the  four  indices  are 
all  equal,  the  quadruple  sum  of  Eq.  (A3)  can  be  written  as  a  single  sum  for  the  case  where  the  indices 
are  all  equal  plus  three  double  sums  for  the  three  cases  where  the  indices  are  equal  in  pairs: 


E{\mR(t)\‘ 


u}  = 


£/74(r  -  kT  -  u) 

k 


+  I  I  (-1  )k~mp2(t-kT-  u) p2{t  —  mT  —  u ) 

k  m^k 

+  £  £  P2(t  -  iT  -  u )  p2(t  -  nT  -  u) 

i  n^i 

+  X  £  ~  iT  ~  u )  P2(t  ~  mT  —  u). 

i  m^i 

Changing  some  indices  and  combining  the  last  three  sums  produce  the  more  compact 

E{\mR{t)\*u)=  J^p'U-kT- u)  +  T  V  [2  +  (^  1  )k~m]  p2{t  -  kT  -  u) p2it  -  mT  -  u). 

k  k  m*k 

This  can  now  be  substituted  back  into  Eq.  (A2)  and  the  integration  brought  inside  the  summations  over 
k.  T 

E{\mR  it)  |4}  =  ±  Zjp4(t-  kT- u)du 

1  k  0 
T 

+  7lJ  £[2+  i-\)k~m]p2it-  kT~  u)p2(t~  mT-  u)  du. 

k  0  m^k 


Applying  the  change  of  variables  x  =  t  -  kT  -  u  in  both  integrals  gives 

f—  kT 

£{|m*(r)|4}=  \  £  f  pHr)dr 

1  k  t—  (jfc+1)  T 
-  t~kT 

+  ^Z  f  T:i2  +  (-l)k-m]p2(T)p2(T  +  (k-m)T)dT. 

1  k  t-(k+\)T  m*k 

It  is  now  clear  that  the  sums  over  k  can  be  eliminated  in  favor  of  infinite  limits  on  the  integrals. 
Changing  the  index  in  the  remaining  sum  to  n  =  k  -  m  simplifies  the  equation  further: 

£(K(f)|4}=  -j  f  p4ir)dr  +  j;  f  £[2  +  (^1)'V(t)p2(t  +  nT)dr. 

— 00  ■*  —00  n^O 

Finally,  the  expression  can  be  compacted  by  combining  the  integrals  into  one  and  expressing  the  first 
term  as  the  missing  n  =  0  term  of  the  sum,  using  a  Kronecker  delta  function  to  get  the  coefficient  right. 

E[\mR  (r)  I4}  =  —  f  £  [2(1—  8„)  4-  (—  1)*]  p2ir)  p2ir  +  nT)dr,  (A4) 

—  00  n 

Equation  (A4)  is  the  general  expression  for  the  fourth  moment  of  the  amplitude  of  this  OQPSK  pro- 
cess  for  real  baseband  pulse  functions. 


Special  Case:  Baseband  Pulse  Functions  With  Limited  Support 

Now,  consider  the  special  case  of  baseband  pulse  functions  (such  as  those  of  Eq.  (18))  satisfying 
the  finite  support  restriction 

pit)  =  0;  \t-  T\>  T.  (A5) 

For  this  class  of  pit)  the  product  p2ir)  p2ir  +  nT)  is  zero  for  all  n  >  2.  The  sum  in  Eq.  (A4)  there¬ 
fore  has  at  most  three  nonzero  terms.  Writing  these  out  explicitly,  distributing  the  integration  over  the 
sum,  and  shrinking  the  limits  of  integration  to  the  smallest  regions  that  might  be  nonzero  produces 
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JT  2  T 

pAW)  dT  +  J  p2(t)p2(t  +  T)  dr  +  J  p2ir)  p2ir  -  T)  dr. 

-  o  T  o  1  T 

A  change  of  variable,  t  =  r-  T,  in  the  integral  third  integral  makes  it  identical  to  the  second.  Carry¬ 
ing  out  this  change  and  splitting  the  first  integral  into  two,  one  from  0  to  T  and  one  from  T  to  27\ 
gives 

T  IT  T 

E{\mR{t)\*)=^\  p*(r)  dr  +  ±f  pH  r)  dr  +  ■ fj  p2(r)p2(r  +  T)  dr. 

1  0  1  T  1  0 

Changing  variables  in  the  second  integral  by  letting  t  =  r  -  T  gives  all  three  integrals  the  same  limits 
and  allows  them  to  be  combined  into  a  single  integral  with  the  integrand  a  perfect  square; 

£{|m*(/)|4}  =  IpHr)  +  pHr  +  T)]2  dr,  (A6) 

1  o 

which  is  the  desired  result  for  pit)  satisfying  Eq.  (A5). 


More  Special  Case:  Baseband  Pulse  Functions  as  Sums  of  Raised  Cosines 


I  now  evaluate  the  fourth  moment  given  by  Eq.  (A6)  for  baseband  pulse  functions  described  by 
Eq.  (18).  I  begin  by  using  Eq.  (18)  to  build  up  the  integrand  of  Eq.  (A6),  simplifying  where  possible. 
From  Eq.  (18)  the  square  of  the  baseband  pulse  function  can  be  written 


P2it)  =  g(t)  £  £2 2akat 

jt-i  /-i 


1  —  cos 


krrt 


1  -  cos 


iirt 


Also  needed  in  the  integrand  of  Eq.  (A6)  is  this  same  quantity  advanced  by  time  T: 


p2it  +  T)  =  git  +  T)  £  £2 2akai 

Ar— 1  /-I 


1  —  COS 


kirit  +  T) 


1  -  cos 


in  it  +  T ) 


Since  cosi<f>  +  kn)  =  (-l)fccos(0),  this  can  be  rewritten 


p2it  +  T)- git  +  T)  £  £2 2akai 

k-\  /-I 


1  -  (— 1)*  COS 


krrt 


1  -  (-1)'  cos 


This  can  now  be  summed  with  Eq.  (A7)  and  the  raised-cosine  products  multiplied  out  to  give 


p2it)  +  p2it  +  T)  =  J  J  2 2aka( 


k- 1  i-l 


2-  [1  +  (-1)*] 


cos 


kirt 


T  -[1  +  (-!)']  cos 


+  [1  +  (— l)fc+']cos  •^1Lt  cos  ~r 


(A7) 


(A8) 


where  the  product  of  the  gate  functions,  not  shown  here,  limits  the  validity  of  this  equation  to  the 
range  0  ^  t  <  T.  Because  this  interval  is  compatible  with  the  limits  of  integration  in  Eq.  (A6),  the 
gating  functions  will  not  be  shown  explicitly  in  the  remainder  of  the  derivation. 

Notice  that  the  second  term  in  the  large  summand  factor  in  Eq.  (A8)  differs  from  the  third  term 
only  in  the  index.  The  total  value  of  the  sum  will  not  change  if  the  index  i  in  the  third  term  is 
replaced  by  k,  yielding  a  slightly  simpler 
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p2(t)+p2(t  +  T)  =  £  £22a*a, 

/-I 


2  —  2  Cl  +  (-l)*]cos  ^ 

+  [l+(-l)fc+/]  cos  cos  — 


(A9) 


The  expressions  of  the  form  [1  +  (-1)*]  that  appear  in  Eq.  (A9)  are  nonzero  only  for  k  even 
This  suggests  the  following  definition: 


even  ( k )  = 


1 ;  k  even 
0 ;  k  odd. 


Using  the  even  function  just  defined,  squaring  Eq.  (A9)  and  substituting  it  into  Eq.  (A6)  results, 
after  some  permutation  of  the  indices,  in 


0  /c— 1  /—l  m— 1  r—  1 


—  28  even(k)  cos 


knt 


+  28  even(k)  even(i)  cos  cos 
+  2 7  evenik  +  /)  cos  cos 

—  28  evenik )  evenii  +  m)  cos  ~r  cos  co§  JHEl 


+  26  even  (£+/)  even  (w+r)  cos  ~~  cos  cos  ^~r  cos 


dt. 


(A10) 


The  products  of  cosines  in  this  equation  can  be  expressed  as  sums  of  cosines  using  the  trigonometric 
identities 

cosa  COS/3  =  2_1tcos(a  +  /3)  +  cos(a  —  /3)  ] , 

cosa  cos/3  cosy  =  2  2[cos(a  +  /3  +  y)  +  cos(a  +  /3  —  y)  +  cos(a  —  (3  +  y)  +  cos(a  —  /3  —  y)  ] , 
and 

cos  a  cos  /3  cos  y  cos  </>  =  2~ 3  [cos  (a  +  /3  +  y  +  </>)  +  cos  (a  +  /3  +  y  -  </>) 

+  cos(a  +  /3  —  y  +  </>)  +  cos(a  +  /3  —  y  —  </>) 

+  cos(a  —  /3  +  y  +  </>)  +  cos(a  —  f3  +  y  —  <f>) 

+  cos  (a  —  /3  —  y  +  <p)  +  cos(a  —  /3  —  y  —  0)  ] . 

Using  the  additional  fact  that 


If  kirt  ,  . 

—  J  cos  — —  dt  =  8k, 
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Equation  (A10)  can  be  integrated  term  by  term  to  yield 


£{|7%(f)|4}  =  XII  Y,akaiamar 

k—  1  /-I  m-1  r— 1 


26 


—  28  even  (k)  8k 

+  2 1  even  {k)  even(i)  (8fc  +  ,  +  §*;_,■) 

+  26  even  (/:  +  /)  (8*+/  +8*;_() 

-  2 6  even  (k)  even  O' +  m)  (8k+i+m  +8fc+/_m  +  8*_,+m  +8fc_,_m) 
+  2 *  even  (k  +  /)  evenim  +  r)  (  8*+,+m+,  +  8/t+H-m-r 

■b  8(t+/-m+r  +  8fc+/_m_f 

“b  8/t—  /+m+r  ~b  ® k—i+m—r 
"b  8 k-i-m+r  +  8 k-t-m-r )  • 


Since  the  indices  are  always  positive,  some  of  the  delta  functions  have  subscripts  that,  because  they 
consist  only  of  positive  terms,  cannot  be  zero.  Eliminating  terms  containing  those  delta  functions  and 
combining  the  two  terms  containing  8jt_,-  gives: 


£’{|w«(t)|4}  =  XII  £  akalamar 
fc-1  /—I  m— 1  r— 1 


26  +  26  [  2  even  (fc)  even  (/)  +  even  (fc  +  /)  ]  8a._,- 


—  26 even(k)  even(/-bm)  (8k+i-m  +  8*_,+m  +  8*_,_m) 
+  23  even  (&  +  /')  even  (m  +  r)  (  8fc+,+m_r  +8k+Hm+r 

+  8fc+/-m-r  +  8*._,+m+,. 
"b  ^ k- i+ m- r~^  ^ k—  i—  m+ r 


+  8 


k—i—m—r 


) 


(All) 


This  can  be  further  simplified  by  removing  some  redundant  references  to  the  even  function.  For  exam¬ 
ple,  consider  the  first  delta  function  in  the  above  equation.  It  is  nonzero  only  if  Ac  —  /  is  zero.  Using 
symbol  implies  for  logical  implication  (if  a  then  b ): 

k  —  i  —  0  =-&■  k  —  i  even 

=■>£►  ( k  odd  A  i  odd )  V  ( k  even  A  i  even ) 

=-^-  -I-  i  even  . 


Therefore, 

even  (fc  +  /)  8*  =  8k  _  ( . 

Similarly, 

even  A  k  ±  i  ±  m  =  0  =—i>  i  +  m  even, 

which,  in  turn,  implies 

even(k)  evenii  +  m)  ( 8k+Hm  +  8k-i+m  +  8/t_,_m)  =  even(k)  (8k+Hm  +  8k-i+m  +  8fc_,_m) . 

The  reference  to  evenim  +  r )  in  Eq.  (All)  can  also  be  shown  to  be  redundant.  Applying  these  sim¬ 
plifications  to  Eq.  (All)  gives 
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E{\mRit)\A}-  ]T  £  £  Z  akaiamar  I  26  +  26  [2  even  (A:)  even(i)  +  1  ]  8k-, 

k- 1  1-1  m-l  r«l  l 

-  26  ( k )  (8k+i.m+8k./+m+8k.,_J 

+  2 3  even  (k  +  0  (  8kH+m_r+8k+l.m+r 

8 k+i—  m-r~^~  8k- t+m+r 

8  k—j+m-r~\~  8 


(A  12) 


The  Mean  of  the  Squared  Amplitude 

I  now  calculate  the  expected  value  of  the  square  of  the  amplitude  of  the  random  OQPSK  process 
described  by  Eq.  (19),  first  with  no  restrictions  on  the  real  baseband  pulse  function  pit)  and  then  with 
the  pit)  of  Eq.  (18)  considered  as  a  special  case  of  the  general  result. 

The  General  Result:  Real  Baseband  Pulse  Functions 


Separate  the  consideration  of  the  random  variable  u  by  formulating  the  desired  moment  in  terms 
of  a  conditional  expectation' 


it)  |2)  =  ~  J  E[\mR  it)  p  u}  du  . 


(A13) 


Now,  calculate  that  conditional  expectation  by  substituting  for  mR  it)  from  Eq.  (19)  and  factoring  all 
the  nonrandom  variables  (including  u,  in  this  context)  out  of  the  expectation: 

E{\mRit)V\u)  =  £  ^E{/3kPm]  jkj~mpit-  kT-  u)  pit  -  mT  -  u) . 

k  m 


The  baseband  pulse  function  pit)  has  been  assumed  real.  Consider  the  remaining  expectation 
EiPkPmT  Because  the  expectation  of  a  product  of  independent  random  variables  is  equal  to  the  pro¬ 
duct  of  the  expectations  of  those  variables,  and  because  the  expectations  of  each  of  the  {/3A}  are  zero, 

E{pk/3m}  =  8m_*. 


Therefore, 


£{|w*(r)|2  u}=  5>2(r-  kT-  u). 


This  can  now  be  substituted  back  into  Eq.  (A13)  and  the  integration  brought  inside  the  summation 
over  k. 

E  {Iw/jW  I2)  =  ^  Z/ P kT~  u^du- 

1  k  0 

Applying  the  change  of  variable  t  =  t  -  kT  -  u  in  the  integral  gives 

t-kT 

E  { I m/f  it)  |2)  =  4  £  J  p2ir)dT. 

1  k  t-(k+l)  T 

It  is  now  clear  that  the  sum  over  k  can  be  eliminated  in  favor  of  infinite  limits  on  the  integral. 

E  {| mR  it)  |2}  =  4  f  P2ir)dT.  (A14) 
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Equation  (A14)  is  the  general  expression  for  this  OQPSK  second  moment  for  real  baseband  pulse  func¬ 
tions. 

Special  Case:  Baseband  Pulse  Functions  With  Limited  Support 

For  the  special  case  of  baseband  pulse  functions  (such  as  those  of  Eq.  (18))  satisfying  the  Finite 
support  restriction 

/>(/)  =  0;  \t-  T\  >  T, 
the  limits  of  integration  can  be  shrunk  to  give 

E  {\mR(t)\2}  =  f  p2(t)  dr.  (A15) 

1  o 


More  Special  Case:  Baseband  Pulse  Functions  as  Sums  of  Raised  Cosines 


I  now  evaluate  the  second  moment  given  by  Eq.  (A15)  for  baseband  pulse  functions  described  by 
Eq.  (18).  Substituting  Eq.  (18)  into  Eq.  (A15)  gives 


E{\mR(t)  |2) 


1  2J  n  « 

tJ  II 


1  —  COS 


knt 


1  cos 


irrt 


dt. 


o  *-i  /-l 

Multiplying  the  raised-cosine  products  and  converting  the  product  of  cosines  to  a  sum  of  cosines  give 


E{\mR{t)V)  =  \,]  £  £«*«, 

1  0  k- 1  /- 1 


a  a  krrt  A  iirt 
4  -  4  cos  — —  4  cos  —jr 


,  ^  (k  +  i)irt  ,  ^  (k  —  i)irt 
+  2  cos  - - — - h  2  cos - - - 


dt. 


The  integration  is  over  integral  multiples  of  the  periods  of  the  cosines,  so  all  the  cosine  terms  disappear 
in  the  integration  except  the  last,  which  has  nonzero  value  only  when  k  —  i  =  0: 

£{K(r)P}=  £  £  u,a/[23  +  228,_/].  (A16) 

fc-i  i-i 


While  the  above  form  will  prove  convenient  for  further  manipulation,  a  more  conventional-looking 
result  can  be  obtained  by  splitting  the  sum  into  two  to  eliminate  the  delta  function 

E[\mR(t) |2)  =  8  £  Y,akai  +  4  £ 

k- 1  i-i  fc-i 


and  then  factoring  the  First  sum: 


£{|m/?(r)|2)  -  8 


+  4  £  al 

k-\ 


The  Variance  of  the  Squared  Amplitude 

The  expressions  given  by  Eqs.  (A16)  and  (A12)  for  the  second  and  fourth  moments  of  the  ampli¬ 
tude  of  the  OQPSK  process  described  by  Eqs.  (19)  and  (18)  can  now  be  combined  according  to  Eq. 
(Al)  to  give  the  variance  of  the  squared  amplitude  of  the  process.  Squaring  the  second  moment  given 
by  Eq.  (A16)  results  in 
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n  n  n  n 


U{KMP})2=  x  I  I  E  akaiamar 

k- 1  i-l  m— 1  r-1 


26  +  2  V;  +  248*_,8m_r 


Combining  this  expression  with  Eq.  (A  12)  in  the  right-hand  side  of  Eq.  (Al)  gives 

n  n  n  n 


Vardw^WR-  HI  I  akaiamar 

At—  1  /-I  m- 1  r-1 


27  eve/z  (£)  eve«  O')  8*_,- 


-248*_,8m_f 

2  even{k)  (8^+/— m  "b  fit— /+m  +8^._/_m) 
+  23  eve«  (*  +  /)  (  8fc+/+m_f  +  8fc+(_m+r 

+  fc+i-m-r  +  $fc-/+m+r 

“t" 


which  is  the  desired  result:  the  variance  of  the  squared  amplitude  of  the  OQPSK  process  described  by 
Eqs.  (19)  and  (18). 
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Appendix  B 

THE  AUTOCORRELATION  AND  SPECTRAL  DENSITY  OF  AN  OQPSK  PROCESS 


In  this  appendix  I  calculate  the  autocorrelation  and  spectral  density  of  the  OQPSK  process  given 
by  Eq.  (19).  The  autocorrelation  of  a  complex  random  process  xit)  is  defined  by 

R^-itu  h)  =  £{x(/i)x*(/2)}- 

Applying  this  definition  to  the  OQPSK  random  process  of  Eq.  (19)  and  using  a  conditional  expectation, 

T 

R  .itht2)  =  ~z,  f  E{mRitx)mRitf)\u)  du  .  (Bl) 

mRmR  T  JQ 


Formulating  the  conditional  expectation  from  Eq.  (19)  yields 

E{mRitx)mRit2)\u)  =  £  'ZE[pkf3m)  jkj~mp(ti-  kT  -  u)pit2- mT  -  u) . 

k  m 

But,  it  is  argued  in  Appendix  A  that  E  {/8*./3m)  =  bm-k.  Therefore, 


E  {mR  it\)mR  it2)  I  u)  =  —  kT  —  u)  p(t2—  kT  u) . 


Substituting  this  back  into  Eq.  (Bl)  produces 


R 


i 

.(r„t2)  =  ^  Xj  p(h-  kT  -  u)p(t2  -  kT  -  u)du. 


mRm*  *  T'o 


A  change  of  variable  v  =  t\  —  kT  —  u  in  the  integral  gives 


t-kT 


R 


(tut2)  =  -^52  /  p(v)  p(v  +  t2- tj  dv  , 


k  l—(.k+\)T 


which  makes  it  clear  that  the  sum  can  be  eliminated  by  extending  the  limits  of  the  integral  to  infinity: 


R 


mRmR 


o° 

(h,h)  =  4  /  p(v) p(v  +  t2~  h)  dv  . 


Since  tx  and  t2  appear  only  in  their  difference,  the  process  mR  it)  is  wide-sense  stationary.  For  such 
processes  it  is  conventional  to  write  the  autocorrelation  as  a  function  of  a  single  variable  representing 
the  difference  between  the  two  times.  Letting  t  =  t2-  h  in  the  above  equation  and  changing  the  vari¬ 
able  of  integration  back  to  the  more  conventional  t , 

oo 

i?ffl  ^  f  pit)  pit  +  t)  dt .  (B2) 

mRmR  T 


This  is  the  general  result  for  the  autocorrelation  of  an  OQPSK  process  with  a  real  baseband  pulse  func¬ 
tion. 


The  spectral  density  of  a  random  process  is  defined  as  the  Fourier  transform  of  its  autocorrelation 
function.  To  transform  Eq.  (B2) ,  first  express  it  as  a  convolution: 

Rm  *(t)  =  -!=  ®p(-t )  ]  • 

mRmR  T 
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Denoting  the  Fourier  transform  of  pir)  by  Fijil),  the  transform  of  pi-r)  is  just  Pi-jCl).  Since 
p(r)  is  real,  Pijn)  is  conjugate  symmetric,  i.e.,  P(j(l)  =  P*(-jtl).  The  Fourier  transform  of  pi-r) 
is  therefore  P*(jil).  Since  the  Fourier  transform  of  the  convolution  of  two  functions  is  just  the  pro¬ 
duct  of  the  Fourier  transforms  of  the  two  functions,  the  desired  spectral  density  can  be  written 

sv;(n>  ~  J  P(ja)P’ija)  -  -L 

The  spectral  density  S  -(D)  is  shown  as  a  function  of  a  rather  than  jCl  to  emphasize  that  when 

mRmR 

integrated  to  obtain  power  it  must  be  with  respect  to  a  only. 


For  the  baseband  pulse  functions  of  interest,  those  given  by  Eq.  (18),  the  (bilateral)  Laplace 
transform  is  given  by  Pis)  =  Gis)Fis)  where  Fis)  is  given  by  Eq.  (16)  or  Eq.  (26)  and  where  G(s) 
is  the  (bilateral)  Laplace  transform  of  the  gate  function  git )  defined  in  Eq.  (5): 


Gis)  = 


2T  e~Ts 


sinh  Ts 
Ts 


The  Fourier  transform  of  pir)  is  obtained  by  evaluating  the  Laplace  transform  Pis)  at  s  =  jil. 
forming  the  substitution  using  Eq.  (26)  for  Fis)  gives 


PijCl)  =  2Tye~JnT  -‘n-  %  - 


n1- 

A:  —  1 


f IT 

Zfc 


2 

T  • 


nr 

kir 


Per- 


and,  consequently, 
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The  above  form  can  be  normalized  to  be  independent  of  the  bit  interval  T  by  expressing  it  as  a 
function  of  a  normalized  frequency  variable  ClT  =  (IT  rather  than  a  function  of  D .  This  must  be 
done  in  such  a  way  that  the  power  in  a  given  frequency  interval,  obtained  by  integration  of  the  density, 
remains  unchanged: 


n,r 


power  (a  i,  a  2)  =  f  s  .(a)  dn  =  f  s  .(a  T)  din  T) 

tn  p  m  p  ^  R  1  1 


Hence, 


5  =  ;S  -(a), 
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or, 
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